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The nature of the nuclear pairing condensate is an active topic of investigation, especially as 
regards its neutron-proton versus identical-particle character, which manifests as the difference be¬ 
tween spin-singlet and spin-triplet pairing. In this work, we probe the recently proposed mixed-spin 
pairing condensates, using a phenomenological Hamiltonian and Hartree-Fock-Bogoliubov theory 
along with the gradient method. In addition to improving the solution of the many-body problem, 
we have calculated a series of physical quantities and examined the robustness of the mixed-spin 
pairing state as the input Hamiltonian is modified. Overall, we find that even though the mixed-spin 
correlation energy is suppressed in comparison to earlier work, the new pairing behavior persists. 
We also discuss the possibility of directly probing the mixed-spin pairing phase. 

PACS numbers: 21.10.-k, 21.30.Fe, 21.60.Jz, 27.60.-|-j 


I. INTRODUCTION 


The interaction between two neutrons is nearly strong 
enough to produce a bound state, a situation that leads 
to fruitful interplay between the physics of neutron-rich 
matter and that of ultracold atomic gases. [T] “Nearly” 
means that while the neutron-neutron scattering length 
is typically much larger than the average interneutron 
spacing, it is not infinite: there is no bound dineutron. 
Similarly, there is no bound diproton. This, added to 
the existence of the deuteron (a bound state of a neutron 
and a proton) led to the proposal of isospin-dependent 
interactions and the conclusion that the isospin-singlet, 
spin-triplet (T = 0, S' = 1) interaction is stronger than 
the isospin-triplet, spin-singlet (T = 1,S = 0) one. 

While the isoscalar (spin-triplet) interaction is stronger 
than the isovector (spin-singlet) one, the pairing appear¬ 
ing in observed nuclei is found between identical particles 
(neutron-neutron and proton-proton). The question of 
what changes in the transition from (two-particle) inter¬ 
action to (many-particle) correlations has therefore been 
actively investigated for a while now. A possible expla¬ 
nation is that the strong nuclear spin-orbit term sup¬ 
presses spin-triplet neutron-proton pairing more than it 
does spin-singlet identical-particle pairing. The concept 
of proton-neutron correlations has therefore been studied 
in a variety of settings, ranging from experiment, binding 
energy systematics, shell-model, and mean-field pairing 
calculations for nuclei, [5HM1 the study of (astrophysi- 
cally relevant) nuclear matter, [53 [55] up to the possibil¬ 
ity that proton-neutron mixing in the particle-hole sector 
is important. m A comprehensive and readable review 
on neutron-proton pairing has appeared recently. [55] 

Another reason why neutron-proton pairing might be 
disfavored (in comparison to identical-particle pairing) is 
that most nuclei don’t have the same number of protons 
and neutrons. This can be easily remedied by eliminat¬ 
ing the isospin polarization, i.e. studying N = Z nu¬ 
clei (as is done in several of the references cited above). 
Combining the two possibilities (reduced effects of the 
spin-orbit field and no isospin polarization). Ref. |29j 


used a phenomenological interaction and Hartree-Fock- 
Bogoliubov (HFB) theory to study heavy nuclei (where 
the surface-to-volume ratio is small) with N = Z, finding 
tantalizing hints of the possibility of spin-triplet pairing, 
which were, however, above the proton dripline. Ref. ISO] 
built on that work by examining N > Z nuclei, finding 
spin-triplet-paired nuclei with unequal neutron and pro¬ 
ton numbers, as well as a smooth transition from spin- 
triplet to spin-singlet, which went through a region of 
“mixed-spin pairing”. The latter nuclei extended below 
the proton dripline and were thus, in principle, accessible 
to experiment. (Intriguingly, a mixed pairing state has 
also been investigated in the context of a one-dimensional 
Hubbard model [8T]). 

In the solution of the HFB equations. Refs. [5^ and 
[30] dropped the Hartree-Fock F term, assuming this was 
effectively absorbed into the model Hamiltonian. In the 
present work, after a brief summary of HFB theory (sec¬ 
tion H AI, the gradient method (section H BI, and our in¬ 
put Hamiltonian (section jll D[ ), we return to this problem 
and explicitly include the F field in our many-nucleon cal¬ 
culations (section [ll E[ ). We then explore the new term’s 
effects on the pairing character for different nuclei (sec¬ 
tion HI A), while also probing particle-number fluctua¬ 
tions (section HI B) and the transition from spin-triplet 
to spin-singlet in more detail (section HI CI. We proceed 
to explore possible experimental signatures of the spin- 
triplet and mixed-spin pairing states (section HI DI. We 
conclude by examining the dependence on the one-body 
(section HIE) and two-body (section HIF) parameters 
in our input model Hamiltonian. 


II. FORMALISM 

A. Hartree-Fock-Bogoliubov theory 

In this section we briefly go over the basics of Hartree- 
Fock-Bogoliubov theory [35], mainly in order to estab¬ 
lish the notation that will be used throughout the paper. 
HFB theory starts from the Fock-space representation of 
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the Hamiltonian: 


^ = E 


ijkl 


( 1 ) 


Then, employing a general Bogoliubov transformation 
one introduces quasiparticle operators pj and Pi'. 

4 = E (2) 
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This already introduces the basic U and V variables, 
which often appear in matrix form as below: 


P 

P^ 


W yt 


(3) 


The Hamiltonian can then be re-expressed in the quasi¬ 
particle basis as follows: 

H = H°°+P'<H^^P+^P'<H^°P'<+ ■■■ (4) 

where the superscripts have the obvious meaning of 
counting the numbers of the quasiparticle creation and 
annihilation operators (the ... represents 
and hermitian conjugates). 

HFB theory amounts to using a variational principle 
and HFB wave functions |$) (along with Thouless’ the¬ 
orem) to minimize the expectation value of H, and then 
writing: 


H = H°° + Y,EkPiPk (5) 

k 

where the E], are known as quasiparticle energies and are 
the eigenvalues in the following HFB equations: 

(S)-(S)®- <*» 


h' A 
-A* -h'* 


where Uk, Vu are the columns of the matrices C/, V. The 
minimization takes place subject to certain constraints, 
such as fixed proton and neutron numbers. Thus, h' = 
h — XiQi, where Qi are constraining fields to be dis¬ 
cussed below and Xi are Lagrange multipliers. The h it¬ 
self consists of the single-particle contribution along with 
a self-consistent field F known from Hartree-Fock theory: 
h = e + T. We thus see that the solution of these equa¬ 
tions rests upon the knowledge of the F and A entities: 


Lzj — ^ ^ VjljkPkl 
kl 

~ 2 ^ ^ ‘^ijklk^'kl 
kl 


(7) 


where the pairing field A encapsulates the physics relat¬ 
ing to a superfluid state. The F and A, in their turn, are 


given in terms of the normal density p and the anomalous 
density k: 

p = V*V'^ 

K = V*U^ (8) 

Note that p is Hermitian, k is skew symmetric, and that 
p and K together uniquely determine the wave function 
| 4 >). 

We see from the Hamiltonian in the quasiparticle ba¬ 
sis, Eq. ([^, that the quasiparticle vacuum, |$), has the 
eigenvalue which turns out to have the form: 

(9) 

This contains contributions from: a) the single-particle 
e term, b) the Hartree-Fock F field, and c) the pairing 
field A. The task of HFB theory can be viewed as the 
minimization of this subject to constraints. It is 
straightforward to see how intimately connected to the 
U and y matrices the solution of the entire problem is. 


B. Gradient Method 


Minimizing the subject to constraints can be a 
daunting task, especially when several constraints are in¬ 
volved. The simplest case of including constraints relates 
to calculations which try to constrain the average particle 
number. In the present work, as in Ref. |d()j . we employ a 
number of constrained parameters to fully investigate dif¬ 
ferent types of pairing condensates, and hence require an 
efficient method to solve the HFB problem. Here we give 
a flavor of the principles involved and refer to Refs. [32]- 
|34l for further details. 

The gradient method begins with a trial wavefunction 
defined by the matrices C/q and Hq, and uses Thou¬ 
less’ theorem to define a wavefunction dependent on an 
antisymmetric matrix Z, |$(Z)): 


mZ)) = A exp I ^ Z,,p\p] I |<i>o), (10) 

i<3 


where A is a normalization factor, Z is the antisymmetric 
Thouless matrix, and Pi annihilates |4)o). The energy 
expectation value, to first order in Z, is then found to 
be: 

E{z) = mz)\Hmz)} 

= AOO + E • (11) 

i<3 

Taking a derivative with respect to Z* of the energy ex¬ 
pectation value with this new wavefunction gives a new 
antisymmetric matrix Zi, namely 

Zu, = mZ)\H MZ)) = 


( 12 ) 
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The parameter t] is somewhat arbitrary, and determines 
how large the step is. This Zi matrix is then used to 
transform the Uq and Vq into new matrices Ui and Vi: 

Ui = Uo + Vo*Z* 

Vi = Vo + [/^Z*. (13) 

which in turn define a new wavefunction |$i). This pro¬ 
cess is repeated until self-consistency, when is zero, 
as in Eq. ([^. 

Adding in constraints is relatively simple. We define 
Qi to be the operator corresponding to the ith constraint, 
which we give a value of Qi. By also defining Lagrange 
multipliers Xi , we can include these constraints by adding 
them into the Hamiltonian: 

H^H-Y^X^Q^. (14) 

i 

This results in the Z\ matrix becoming: 


reformulated in the HFB formalism, using the ground 
state 14)). Thus, the variance becomes 

(A6)2 = ^Tr (0^0002^ ^ ('20) 

Below we will use separate operators for the proton and 
neutron numbers and will thereby determine the proton 
and neutron number fluctuations (defined as the square 
root of the variance). 

D. Hamiltonian 

We consider a phenomenological Hamiltonian, consist¬ 
ing, as usual, of a one-body part and a two-body inter¬ 
action term, represented as 

ij ijkl 



In order to determine what the Xi are, we impose the con¬ 
dition that the expectation of Qi equals Qi, to first order 
in Z. The operators also have their own quasiparticle 
basis representations, similar to that of the Hamiltonian, 
with their corresponding Q™, and Qf°, so they have 
an expectation value of 


ft) = QT+Y iKkiQl%i + ZixiQfu) = qT+ZvQT- 


k<l 


Our condition then requires Zi 


(16) 

Qf° to be zero. Using 


our new Zi from Eq. (15) with the included constraints 
then leads to a system of equations involving the Xf. 




This allows us to relatively easily solve for Xi. Once these 
Ai’s are determined, this Zi is used as normal. In the 
present work, up to 8 constraining fields are explored, 2 
for the neutron and proton particle numbers and 6 for 
the pairing configuration amplitudes. 


C. Particle Number Fluctuations 


v^Jkl = j d^rSr'x*{v)x*{v')V{v,v')xk(x)xii^') (21) 

the only difference from Eq. Q being that now we are not 
using antisymmetrized matrix elements. The Xi(r) are 
coordinate-space wave functions associated with a single¬ 
particle Hamiltonian. 

In order to probe the essential features of pairing in 
heavy nuclei, we here follow Refs. jlHl ISO] and take the 
one-body part e to contain a kinetic energy, a Woods- 
Saxon well, and a spin-orbit term: 

s= -h Vwsf[r) + 1 • cr Vso-^ — 

Zm r dr 

= 1 + fXr-R)la (22) 

where R is the radius of the nucleus, R « 1.27A^/^. The 
other parameters are taken to be: a = 0.67 fm, Vws = 
—50 MeV, and Vso = 33 MeV fm^, [29l[30j though these 
are also probed in more detail in a later section. 

Similarly, the two-body interaction term is taken to be 
of contact form: 

6 

V (r, r') = Y (r - r')PL=oPa 

a —1 

= I (^Vt -I- Us -I- {vt - Vs)(T ■ a') 


One can also use the above formalism to examine the 
question of particle-number fluctuations in this theory. 
[34] For a single-particle operator Q: 

\ , 

( 18 ) 

we can calculate the statistical variance 

(aQ2) = (Q2) _ (0)2 (19) 


x<5(3)(r-r')PL=oj (23) 

The operator Pl=q projects onto states with zero total or¬ 
bital angular momentum. This restriction means nuclear 
deformations are not considered, but it simplifies the cal¬ 
culations involved by imposing a block-diagonal struc¬ 
ture on the U and V matrices. The index a represents 
one of the six pairing configurations shown in Table |HD[ 
which are explicitly taken into account when we construct 
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a 

1 2 3 4 5 6 

{S,Sz) 

{T,Tz) 

(0,0) (0,0) (0,0) (1,1) (1,0) (1,-1) 
(1,1) (1,0) (1,-1) (0,0) (0,0) (0,0) 


TABLE I. Pairing channels organized by spin-isospin. S is 
the two-nucleon spin, Sz its a-projection; T is the two-nucleon 
isospin, Tz its ^-projection. 


the pairing condensates. The operator projects onto 
a particular pairing configuration. Obviously, the pair¬ 
ing strengths Va can take on only two values, Vs and Vt 
for spin-singlet and spin-triplet, respectively, as explicitly 
shown in the second step of Eq. (23). These are taken to 
have values of 300 MeV and 450 MeV, respectively, to be¬ 
gin with (having been fit to shell-model matrix elements) 
though their values will be varied in a later section. Fi¬ 
nally, the Dirac delta interaction is smeared by including 
only orbitals with single-particle energies ±5 MeV away 
from the Fermi energy. 


E. Implementation of P and A 


In Refs. [291 [30], given that the Hamiltonian was phe¬ 
nomenological, the r matrix was assumed to be implic¬ 
itly included in order to simplify the calculations. In this 
work, we explicitly introduce the T matrix correspond¬ 
ing to the specific interaction described in the previous 
section. 

If we introduce the interaction of Eq. (23) into the 


expression for the two-body matrix elements in Eq. (21), 
we can express things more compactly: 


'^ijkl — / , 


, J (r-r')(fj|r, r')P„Pi^o(i’, r'\kl). 

(24) 

The kets \ij) contain all spatial, angular momentum, 
spin, and isospin quantum numbers. We can thus sepa¬ 
rate them into a spatial ket (which does not affect) 
and another ket containing the remaining quantum num¬ 
bers (which Pa does affect). Expressing the projection 
operator Pa as |a) (a| allows us to introduce the follow¬ 
ing factors Aa^ij, naturally defined in terms of Clebsch- 
Gordan coefficients: 


= 72 


1 


1 




SaSz 
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2 ) r) I "ZJ 




X (—1)^’ 


(25) 


If we now integrate over r', write the spatial bra-kets 
in terms of single-particle wave functions, extract their 
spherically symmetric part, and carry out the angular 
integral we find: 


1 
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Vijki = / drr^(f)*{r)(f)*Ar)(f)k{r)4'i{r)Aa,jrAa,ki- 


47r 


( 26 ) 


where the (piir) are radial functions and we also switched 
to antisymmetrized matrix elements. 

It is now straightforward to calculate the A and T 
fields, starting from Eq. ([^. For A we have: 

Aij = — ^ ^ 'jJjjkjk^kl 
kl 

6 

= 

Ct — 1 


J drr^(l)*{r)(l)*{r)Aa,ijK.a{r) (27) 


where 


^ X! (28) 

^ k<l 

as shown in the Appendix of Ref. |29j . Similarly, we now 
find for the T matrix: 

~ ^ ^ ^iljkPkl 
kl 

= ^ ^ f drr'^(j)*{r)(j)*{r)(j)k{r)(|)l{r)Aa,^lPlkAa,kj■ 

(29) 

For both the A and T terms, we have written the final re¬ 
sults in a form that highlights the matrix multiplications 
involved. Note that since there is no trace of a matrix 
multiplication (as in the case of the A matrix) there is 
no separation of steps needed for the T matrix. 


III. RESULTS 

A. Pairing below the N = Z line 

In what follows, we will be interested in probing the 
effects of pairing in heavy nuclei. To do that, we calculate 
the ground-state energy of a many-nucleon system when 
no pairing is present (calling that Eq), as well as the 
energy when pairing is turned on (calling that E). To 
ensure that Eq corresponds to the case where no pairing 
is present, we explicitly constrain all 6 pairing amplitudes 
(kq = PiPaiA) to zero. The difference between the two 
energies is the correlation energy. 

Ecorr = Eq - E (30) 

In Fig.[2 nuclei below the N = Z line have been mapped 
out for neutron numbers between 50 and 75. (Here 
and in the rest of this paper results shown include the 
r term in the many-body calculations.) Some nuclei 
have very small correlation energies, indicated by blank 
spaces, while most of the remaining nuclei exhibit spin- 
singlet pairing, denoted by green squares. Between neu¬ 
tron numbers 60 and 70 there is an island of nuclei that 
exhibit spin-triplet pairing (red diamonds), and another 
region of mixed pairing (blue circles). The mixed-spin 
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FIG. 1. (Color online) Chart of nuclides with N > Z from 
A'’ = 50 to = 75. Blank squares indicate nuclei with no 
pairing {Ecorr < 0.5 MeV). Green squares denote a pair¬ 
ing condensate of mostly spin-singlet character, red diamonds 
those of a spin-triplet character, and blue circles indicate a 
mixture of the two types of pairing. 


pairing has been defined to be a state for which the spin- 
singlet amplitude is between one quarter and three quar¬ 
ters of the total pairing amplitude. Also shown in the fig¬ 
ure (black dashed line) is the proton dripline, produced 
based on Ref. [55] . 

Comparing this figure with Fig. 1 in Ref. [SD], it can 
be seen that the inclusion of the F matrix has a slight 
shifting or suppressing effect on spin-triplet pairing, as 
both the spin-triplet and mixed-spin regions have shifted 
closer to the N = Z line. Nevertheless, the majority of 
the mixed-spin region remains below the proton-dripline, 
and thus may still be relevant to experiment. Detailed 
values of the correlation energy in the cases of including 
and excluding the F term will be discussed in connection 
with Fig. below. 


B. Number fluctuations 


Second, there are two areas of relatively higher fluctua¬ 
tions in both numbers, one centered at {N, Z) = (70, 60) 
and a smaller one at {N,Z) = (57,55). As could have 
been expected, these closely correlate with regions having 
large correlation energy. 


C. Singlet to triplet transition 


As in Ref. [5D|, in order to more carefully examine 
the transition from spin-triplet to spin-singlet as a func¬ 
tion of increasing N — Z, we study selected nuclei with 
A = 132. This line was chosen as it exhibits examples of 
the three possible states: spin-singlet, spin-triplet, and 
mixed pairing condensates. Specifically, the three nuclei 
were chosen to be go^Nd, e 4 ^Gd, and gg^Dy. In addition 
to the ground-state results shown in Fig. [l] we have here 
artificially constrained the spin-triplet and spin-singlet 
amplitudes away from the minima (explicitly: channels 
5 and 1/3 from Table IID). The resulting correlation en¬ 
ergies are shown on a contour plot in Fig. [^ 

The plot for gg^Nd on the left has peaks at (±14,0), 
indicating a spin-singlet condensate, which is the more 
ordinary state of affairs. The unconstrained correlation 
energy at the minimum is ~ 6.5 MeV, to be compared 
with the ^ 8 MeV when the F term is excluded. The plot 
for gg^Dy on the right shows peaks at (0, ±22), indicating 
a spin-triplet condensate. Here the correlation energy is 
~ 4 MeV, to be compared with ^11 MeV without the F 
term. This is a major reduction, clearly reflecting the fact 
that the F term serves to suppress the exotic spin-triplet 
state. Finally, the middle plot for H^Gd shows peaks at 
(±8,±19), indicating a mixed-spin pairing condensate. 
In this case the correlation energy is ~ 4 MeV, to be 
compared with ^ 7 MeV without F, an intermediate re¬ 
duction. Overall, as N — Z increases, the type of pairing 
does not suddenly switch from one to the other; rather, 
it gradually changes from spin-triplet to spin-singlet, a 
conclusion that survives from the F-less case. Note, how¬ 
ever, that the correlation energies in all three cases have 
been reduced (though not sufficiently to drown out the 
effects under study). 


Since HFB theory does not conserve particle number 
exactly, but only on average, it is worthwhile to inves¬ 
tigate the particle-number fluctuations. To do this, we 
have employed Eq. (20), with separate operators for the 
proton and neutron numbers. The results are shown in 
Fig. 15 Overall, we see that the fluctuations never exceed 
4 particles (or 3%). 

Two main features stand out in this figure. First, there 
are regions of the chart that have zero particle-number 
fluctuation. As can be seen by comparing with the cor¬ 
relation energies that led to Fig. no fluctuations are 
trivially equivalent to the absence of pairing, which au¬ 
tomatically restores good particle number. There is a 
pretty consistent line at the magic number oi Z = 50. 


D. Pairing gaps 

Up to this point we have examined correlation energies, 
which reflect properties at the many-particle level. It is 
worthwhile to investigate possible single-particle proper¬ 
ties. One way to do so is by calculating the pairing gap 
from the odd-even mass differences. This can be arrived 
at via a simple difference formula: 

(n) = E{n) - i \E{n - 1) + U(n ± 1)], (31) 

where n is an odd neutron or proton number, with the 
other nucleon number set to an even value, and E(n) is 
the ground state energy of the nucleus. 
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FIG. 2. (Color online) Heat maps showing the fluctuations in neutron number (left panel) and proton number (right panel) for 
heavy nuclei on and off the N = Z line. They range from 0 particles to a little over 3 particles. 
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FIG. 3. (Color online) Contour plots showing the correlation energy (in MeV) in 3 different nuclei on the A = 132 line as a 
function of spin-singlet and spin-triplet pairing amplitudes. From left to right are epNd, gl^Cd, and ee^Dy. eo^Nd exhibits 
mostly spin-singlet pairing, and ge^Dy exhibits mainly spin-triplet pairing. gt^Gd exhibits a mixed-state pairing that shows a 
smooth transition from spin-triplet to spin-singlet as N — Z increases. 


The results are shown in Fig. As in Ref. m, we find 
a signature of spin-triplet pairing in a chain of nuclei one 
unit away from N = Z (above the proton drip line). 
What’s different here is that we also find a sequence of 
small pairing gaps on the A = 133 line. This results 
from a more thorough search, imposing the odd-number 
parity per block across blocks, as well as the effects of 
the r matrix. The placement of the A = 133 line is not 
coincidental: it reflects the fact that above that line new 
states become available. Additionally, since orbitals are 
included in the HFB space if their energies lie between 
the bounds E^nin and ifmax, changing the energy bounds 
also changes the energy levels available for calculations 
for each mass number. We have checked the effects of 
increasing the AE and found no qualitative changes. 


E. One-body potential parameters 


Our discussion of pairing gaps led to a study of the 
dependence (or lack thereof) of observed effects on spe¬ 
cific model parameters. Continuing in that vein, we also 
varied the one-body and two-body parameters appearing 
in Eqs. (22) and (23), in order to determine the sensitiv¬ 


ity of the different pairing condensates to details of our 
input Hamiltonian model. 

We start by ad hoc varying the Vso one-body parame¬ 
ter, which controls the strength of the spin-orbit term. 
Qualitatively, based on the ideas discussed in the In¬ 
troduction, we expect that reducing Vso would lead to 
a strengthening of spin-triplet pairing or, reversely, in¬ 
creasing Vso would lead to a suppression of spin-triplet 
pairing. Figure shows the results of carrying out such 
calculations. Since our value for the spin-orbit strength 
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above was Vso = 33 MeV fm^, we here depart by 7 MeV 
fm^ on either side. The nuclei investigated are the same 
as those in Fig. ^ from left to right, they are gg^Nd, 

132pj 132 t^,, 

64 ’-'Qj 66 

The results for gg^Nd on the left show a fairly steady 
increase in correlation energy as Vso increases. The na¬ 
ture of the pairing also remains steadily spin-singlet. Ad¬ 
ditionally, we determined that the inclusion of the T ma¬ 
trix once again causes a slight drop in the correlation 
energy, approximately 25%. The results for If Dy on the 
right are similar, in that the spin-triplet character is not 
violated. Of particular interest is the decrease in corre¬ 
lation energy as Vso increases, rather than an increase. 
This reflects the qualitative expectation that increased 
Vso is not favourable to spin-triplet pairing. The results 
for If Gd, shown in the middle of Fig. are the most 
interesting. Similar to the |f Nd case, the correlation 
energy experiences an increase as Vso increases. But, 
also as Vso increases, the pairing nature changes from 
spin-triplet to a mixed-pairing state. 


From general considerations one expects that the pres¬ 
ence of the spin-orbit held (pushing a partner to higher 
energies) suppresses pairing. This suppression is active in 
both the isovector and isoscalar channels, see e.g. Fig. 11 
in Ref. |28]. On the other hand, the left panel of Fig. 
shows the correlation energy increase as the spin-orbit 
strength is increased. We have traced this to the orbitals 
selected in the A=132 calculations, along with the rele¬ 
vant mixing. It’s worth noting that the trend exhibited 
by the correlation energy in |f Nd is analogous to the 
behavior of the single-particle contribution to the total 
energy: the spin-orbit field has a more pronounced effect 
here than in the mixed-spin or spin-triplet nuclei. 

Overall, Fig. [^supports the proposed suppressing na¬ 
ture of the nuclear spin-orbit interaction on spin-triplet 
pairing. Importantly, it shows that the spin-triplet and 


N 



FIG. 4. (Color online) Pairing gaps for 55 < A < 74 and N > 
Z. Most nuclei have intermediate-size pairing gaps (green 
squares), though two lines exhibit reduced gaps (blue circles). 


Nucleus 

Vws 

Vso 

-^corr 

singlet fraction 

IfNd 

-50 

33 

6.432 

1 


-50 

33.28 

6.482 

1 


-47 

33.28 

8.016 

1 

M"Gd 

-50 

33 

3.904 

0.2784 


-50 

34.75 

4.209 

0.3027 


-49 

34.75 

4.138 

0.3096 

efDy 

-50 

33 

4.067 

0 


-50 

35.48 

3.894 

0 


TABLE II. Correlation energies for |f Nd, |f Gd, and |f Dy 
with various Vws and Vso. Also shown is the “singlet frac¬ 
tion” , namely the ratio of the spin-singlet amplitude and the 
total pairing amplitude. 


mixed-spin states, while depending on the overall magni¬ 
tude of Vso 1 do not result from a “magic” choice that is 
finely tuned: both exotic states appear for a large spec¬ 
trum of possible spin-orbit strength values. 

The logical extension of the trends shown in Fig. 
is to examine what would happen if in our calculations 
we switched off the spin-orbit strength, Vso j completely. 
This is, of course, unphysical, so functions only as a con¬ 
sistency check. The results are shown in Fig. we find 
clear signatures of spin-triplet pairing across the entire 
N = Z line, as could be expected from the qualitative 
arguments given in the Introduction, as well as slightly 
off it. This figure also nicely illustrates the other cause 
behind the suppression of spin-triplet pairing: the isospin 
asymmetry gradually turns spin-triplet states into mixed- 
spin states (further below the N = Z line) which even¬ 
tually turn into spin-singlet states. 

In the spirit of more carefully probing the input Hamil¬ 
tonian parameters, we have also considered the possibil¬ 
ity that these may change according to how neutron-rich 
the nucleus under study is. In Ref. [30] the parame¬ 
ters Vws and Vso were kept constant regardless of the 
value oi N — Z or A. Here, motivated by a standard 
formula, [SB] we also examine the effect of changing these 
parameters (individually or in concert) on the correlation 
energies for Ig^Nd, H^Gd, and ||^Dy. In keeping with our 
default values, we explore the following dependences: 


VWs = ( -50 -H 33 ^^ MeV, 
Vso =rl ( 22 - MeV, 


(32) 


where rg = 1.27 fm. The results are shown in TablejTTj 
First, we note that the singlet fraction (defined as 
the ratio of spin-singlet and total pairing amplitudes) 
is not impacted by these alternative choices for Vws ^nd 
Vso: remains firmly on the spin-singlet side, H^Gd 

still exhibits mixed-spin pairing, and |g^Dy is still over¬ 
whelmingly a spin-triplet nucleus (of course, the main 
point of Eq. (321 was to explore an N — Z dependence. 


which would not show up for |g Dy). Next, we notice 
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FIG. 5. (Color online) Correlation energies (in MeV) vs. varying Vso (in MeV fm^). These are plotted for the same nuclei as in 
Fig. i epNd in the left panel, et^Gd in the middle panel, and el^Dy in the right panel. As in Fig.[^ red indicates spin-triplet, 
green spin-singlet, and blue a mixed-state. 


that changing the Vso parameter alone by these small 
amounts does not have an appreciable effect on the cor¬ 
relation energy or on the nature of the pairing conden¬ 
sate: this is hardly surprising, since precisely the same 
modification was carried out for the purposes of Fig. 
Moving on to the lines of the table showing a change of 
both Vws and Vso- shows an increased correla¬ 

tion energy, though no fundamental change in the nature 
of the pairing, while gt^Gd is barely impacted. 


F. Pairing strength variation 


As noted in section IID 


Ref. 


the pairing 


strengths Vs and Vt (for spin-singlet and spin-triplet, 
respectively) were arrived at by comparing with phe¬ 
nomenological shell-model Hamiltonians. The Vg was 
estimated to be ~ 280 and the ratio between the two 
parameters was found to be ~ 1.65, so the two parame- 



50 55 60 65 70 75 


N 

FIG. 6. (Color online) Chart of nuclides with N > Z from 
V = 50 to V = 75 for Vso ~ 0. Notation is as in Fig. A 
vanishing spin-orbit strength leads to the appearance of many 
new spin-triplet and mixed-spin nuclei. 


ter values were conservatively estimated to be 300 MeV 
and 450 MeV, respectively. These were also the strengths 
used in Ref. |30j and in the present work up to this point. 
Here, once again, we attempt to see if there is a sensitive 
dependence on the specific values chosen or, alternatively, 
if the conclusions on the pairing condensate character(s) 
are fairly robust. We have examined Vg values from 250 
to 350 MeV, for vt/vg ratios of 1.25, 1.4, 1.5, 1.6, and 
1.75. The results for the correlation energies are shown 
in Fig. As in earlier figures, on the left is on 

the right is gg^Dy, and in the middle is ifGd. 

An unmistakeable trend is that as we move to the right 
of the x-axis in each figure, increasing both Vg and Vt, the 
correlation energy also increases, a natural result of hav¬ 
ing stronger interactions. Note that for gg^Nd the spread 
of the results for different Vt/vg ratios is relatively small: 
for most strengths and ratios the correlation energy is ap¬ 
proximately the same. Even for this clearly spin-singlet 
paired nucleus, a very large Vt/vg ratio starts to lead to 
mixed-spin pairing. The behavior of If Gd and if Dy is 
quite different: for these nuclei different Vt/vg ratios lead 
to considerably different correlation energies. For small 
ratios both nuclei exhibit spin-singlet pairing, while over¬ 
all If Gd spans the gamut from spin-singlet, to mixed- 
spin, to spin-triplet pairing: the presence of mixed-spin 
pairing for reasonable values of the pairing strengths is 
fairly well established here. In the case of |f Dy, we see 
that for small pairing strength ratios (at most values of 
Vg) there is no pairing whatsoever (zero correlation en¬ 
ergy), while for any other ratio we find a clean signal 
corresponding to spin-triplet. 


IV. SUMMARY AND CONCLUSION 


In summary, we have improved the solution of the 
Hartree-Fock-Bogoliubov problem in comparison to ear¬ 
lier work, by employing the gradient method and explic¬ 
itly including the Hartree-Fock F field. We have then car¬ 
ried out calculations for nuclides with N > Z ior N = 50 
to N = 75. These were carried out by constraining the 
average proton and neutron particle numbers. We find 
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FIG. 7. (Color online) Correlation energies (in MeV) vs. varying Va, with Vt/Vs ratios of 1.25, 1.4, 1.5, 1.6, and 1.75. On the 
left is go^Nd, on the right is gg^Dy, and in the middle is g|^Gd. All plots indicate a change from spin-singlet (green) towards 
spin-triplet (red) possibly with an intermediate mixed-state (blue) as the Vt/Va ratio increases, and an increase in correlation 
energy as Vs increases. 


that including F leads to a suppression of the correla¬ 
tion energy, which is not, however, sufficient to elimi¬ 
nate the mixed-spin pairing condensate for a number of 
nuclei below the proton-drip line. We have also taken 
the opportunity to examine the particle number fluctu¬ 
ations, which end up being reasonably small and closely 
following the magnitude of the correlation energy. Fur¬ 
thermore, we have investigated the effect of artificially 
constraining the spin-singlet and spin-triplet pairing am¬ 
plitudes away from the ground-state minima: the results 
are qualitatively unchanged. Regarding the pairing gaps, 
we found a sequence of spin-triplet character (above the 
drip line) as well as a line of reduced pairing gaps below 
the proton drip line. We have also attempted to modify 
our model Hamiltonian one-body and two-body parame¬ 
ters. We found that the spin-orbit strength quenches the 
spin-triplet pairing, though not dramatically so for values 
close to our default ones. Similarly, the pairing strengths 
impact the nature of the pairing condensate, but we see 
no signs of our default values being finely tuned. 

In terms of future work, the natural next step would 
be to extend our mean-field theory so that it can han¬ 
dle broken symmetries, like particle number and angular 
momentum. m For example, at this stage we have not 


made statements on the ground-state spins for the nu¬ 
clei that we find exhibit mixed-spin pairing. Such an 
extension would allow one to make predictions for spec¬ 
troscopic quantities as well as two-particle transfer reac¬ 
tions. Another natural avenue of future work would be to 
address nuclear deformation [24) : while the prediction of 
mixed-spin pairing for experimentally accessible nuclei is 
tantalizing, the results we have produced have been until 
now restricted to the case of spherical symmetry, which is 
known not to apply to the region of interest. The pairing 
would certainly be weakened in a fully three-dimensional 
Hartree-Fock-Bogoliubov calculation, though it remains 
to be seen by how much and in what way. 
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